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ABSTRACT 

Context. Gravitational collapse theory and numerical simulations suggest that the velocity field within large-scale galaxy filaments is 
dominated by motions along the filaments. 

Aims. Our aim is to check whether observational data reveal any preferred orientation of galaxy pairs with respect to the underlying 
filaments as a result of the expectedly anisotropic velocity field. 

Methods. We use galaxy pairs and galaxy filaments identified from the Sloan Digital Sky Survey data. For filament extraction, we 
use the Bisous model that is based the marked point process technique. During the filament detection, we use the centre point of each 
pair instead of the positions of galaxies to avoid a built-in influence of pair orientation on the filament construction. For pairs lying 
within filaments (3012 cases), we calculate the angle between the line connecting galaxies of each pair and their host filament. To 
avoid redshift-space distortions, the angle is measured in the plain of the sky. 

Results. The alignment analysis shows that the orientation of galaxy pairs correlates strongly with their host filaments. The alignment 
signal is stronger for loose pairs, with at least 25% excess of aligned pairs compared to a random distribution. The alignment of galaxy 
pairs and filaments measured from the observational data is in good concordance with the alignment in the Millennium simulation 
and thus provides support to the ACDM formalism. 
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1. Introduction 

The main driving force of structure formation in the expand¬ 
ing Universe is the gravitational collapse. The collapse ampli¬ 
fies any initial anisotropy of matter distribution. As first pointed 
out by Zel’dovich (1970), the contraction is locally dominated 
by one dimension at a time. The speed of the collapse depends 
on the local density and the initial amplitude of perturbations, 
thus in numerical simulations, a complex pattern of nearly 2- 
dimensional sheets, nearly 1-dimensional filaments and roughly 
spherical small density peaks (knots) develops, framing large 
under-dense void regions (Doroshkevich et al. 1980; Klypin & 
Shandarin 1983; Davis et al. 1985; Gramann 1993; Sheth 8c van 
de Weygaert 2004, etc). Within this network, matter flows along 
well defined paths, with the voids feeding the sheets, the sheets 
feeding the filaments and the filaments finally channeling ev¬ 
erything towards the knots (Cautun et al. 2014, and references 
therein). 

Observational confirmation for the existence of a characteris¬ 
tic pattern of large-scale inhomogeneities came already decades 
ago (Joeveer & Einasto 1978; Joeveer et al. 1978; Peebles 1980; 
de Lapparent et al. 1986; Bond et al. 1996, etc). However, the 
low number density of sufficiently luminous galaxies sets severe 
limits to the possibilities of comparing the actual structures to the 
simulated ones. As a result, the fine structure and the dynamics 
of the actual cosmic web can be directly probed only within our 
local volume, which may not be representative of the Universe at 
larger scales. Nevertheless, the existing peaces of evidence about 
the dynamics within our super-galactic neighbourhood support 
the void-sheet-filament-cluster sequence of matter flows (Tully 
et al. 2014; Tempel 2014; Karachentsev et al. 2014). 

As an alternative to the simple mapping of structures, a 
method for testing and constraining structure formation scenar¬ 


ios is provided by the alignment of structures. Anisotropic col¬ 
lapse and tidal torques should induce spinning up and stretching 
out of dark matter haloes along the general velocity field within 
a large-scale structure element (Codis et al. 2012; Aragon- 
Calvo & Yang 2014; Libeskind et al. 2014; Welker et al. 2014). 
Observationally, corresponding hints for some alignment be¬ 
tween various structures have been found: elliptical galaxies and 
satellite systems around massive galaxies tend to be elongated 
in the direction of galactic filaments (Tempel et al. 2013; Lee 
& Choi 2015; Zhang et al. 2015), while galaxy discs tend to be 
oriented perpendicular to filaments (Tempel et al. 2013; Tempel 
& Libeskind 2013). The former findings can be explained with 
galaxy mergers and satellite accretion occurring along the fila¬ 
ments, while the latter effect can be reconciled with matter infall 
from voids and sheets around the filaments (Dubois et al. 2014). 


In this paper we take a new step in the alignment probing by 
measuring the angle between the orientation of galaxy pairs and 
the orientation of their host filaments. The pairs and filaments 
have been extracted from the Sloan Digital Sky Survey (SDSS). 
Eor filament detection we use the Bisous model, which is based 
on a marked point process (Tempel et al. 2014a). The Bisous 
filaments have already proven to be handy for probing the galaxy 
content and other properties of the cosmic web (e.g. Tempel et al. 
2014b,c; Guo et al. 2015). 


Throughout the paper we assume the Wilkinson Microwave 
Anisotropy Probe (WMAP) cosmology: the Hubble constant 
Hq - 100 kms^'Mpc^', with h — 0.697, the matter density 
Dm = 0.27 and the dark energy density Qa = 0.73 (Komatsu 
etal. 2011). 
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Fig. 1. The distribution of pair separations in the plane of the 
sky. The vertical dashed line splits the sample into two equal 
subsamples. 

2. Data and methods 

2.1. SDSS galaxy and group samples 

This work is based on the SDSS data release 10 (York et al. 2000; 
Ahn et al. 2014). We use the galaxy and group samples as com¬ 
piled in Tempel et al. (2014d), which cover the main contigu¬ 
ous area of the survey (the Legacy Survey). Galaxy groups are 
extracted using the conventional friend-of-friend (FoF) grouping 
algorithm with a varying linking length (see Tempel et al. 2014d, 
for more details). The flux-limited catalogue (apparent r-band 
magnitude limit is 17.77) extends to the redshift 0.2 (820 Mpc) 
and includes 588193 galaxies and 82458 groups with two or 
more members. 

In Tempel et al. (2014d) the redshift-space distortions of 
galaxy groups are suppressed using the rms sizes of groups in 
the plane of the sky and their rms radial velocities as described in 
Liivamagi et al. (2012). Galaxies in groups are ascribed new ra¬ 
dial distances in order to make the spatial distribution of galaxy 
groups approximately spherical. Such a compression removes 
the artificial line-of-sight filament-like structures and was neces¬ 
sary to retrieve a realistic filamentary pattern in the cosmic web. 

For our analysis we use only galaxy pairs (galaxy groups 
with two members in the Tempel et al. (2014d) catalogue) that 
are located within filaments (see Sect. 2.3 for the filament cat¬ 
alogue description). We consider a galaxy pair to be located in 
a filament if the geometrical centre point of the pair is closer 
than 0.5 Mpc to the filament axis. In total, this limits our sample 
to 3012 pairs. To avoid redshift space effects, the separation be¬ 
tween the galaxies in a galaxy pair is measured in the plane of 
the sky. The distribution of pair separations for pairs located in 
filaments is shown in Fig. 1 . 

Due to fibre collision in SDSS, some of the pairs in our sam¬ 
ple are actually groups with three or more members. According 
to Tempel et al. (2012) the fraction of pairs that are affected by 
fibre collision is about 5%, hence, we can ignore this effect in 
the current analysis. 

2.2. Seml-analytic galaxies In simulation 

For a quick interpretation of the results, we use the distribu¬ 
tion of galaxies in the Millennium simulation (Springel et al. 
2005), in which dark matter haloes are populated with galaxies 
according to a semi-analytical galaxy formation model galform 
(Bower et al. 2006). The simulated galaxy data and the galactic 
filaments are the same as used in Guo et al. (2015) and Tempel 
et al. (2015) - we refer to these papers for a more detailed de¬ 
scription. For the analysis, we only have considered galaxies 
brighter than Mr = -19. The galaxy pairs (groups with two 
members) were extracted using the FoF algorithm with a link¬ 
ing length 0.5 /; 'Mpc. The pairs are separated into equal-size 


close and loose pair samples according to the real separation be¬ 
tween galaxies. The limiting distance is 0.2 /z^'Mpc (i.e. roughly 
the same as for SDSS). 


2.3. Galaxy filament sample for the SDSS data 

The detection of filaments is performed by applying an ob¬ 
ject/marked point process with interactions (the Bisous process; 
Stoica et al. 2005) to the spatial distribution of galaxies. This al¬ 
gorithm provides a quantitative classification that complies with 
the visual impression of the cosmic web and is based on a robust 
and well-defined mathematical scheme. More detailed descrip¬ 
tions of the Bisous model can be found in Stoica et al. (2007, 
2010) and Tempel et al. (2014a). For reader convenience, a brief 
and intuitive summary is provided below. 

Our model approximates the filamentary web by a random 
configuration of small segments. We assume that locally, galaxy 
conglomerations can be probed with relatively small cylinders, 
which can be combined to trace a filament if the neighbouring 
cylinders are oriented similarly. An advantage of the approach is 
that it relies directly on the positions of galaxies and does not re¬ 
quire any additional smoothing for creating a continuous density 
field. 

The solution provided by our model is stochastic. Therefore, 
we find some variation in the detected patterns for different 
Markov chain Monte Carlo (MCMC) runs of the model. On the 
other hand, thanks to the stochastic nature of the method we gain 
a morphological and a statistical characterisation of the filamen¬ 
tary pattern simultaneously. 

In practice, after fixing an approximate scale of the filaments, 
the algorithm returns the filament detection probability field to¬ 
gether with the filament orientation field. Based on these data, 
filament spines are extracted and a filament catalogue is built in 
which every filament is represented by its spine: a set of points 
that define the axis of the filament. 

Our spine detection follows two ideas. First, filament spines 
are located at the highest density regions outlined by the filament 
probability maps. Second, in these regions of high probability 
for the filamentary network, the spines are oriented along the 
orientation field of the filamentary network. See Tempel et al. 
(2014a) for more details of the procedure. 

For the current study, we have rebuilt the filament catalogue 
in order to eliminate the potential influence of galaxy pairs to the 
filament orientation. For filament extraction, each galaxy pair is 
replaced with a single point (the centre point of the two galaxies 
in the 3D space). Doing so, the orientation of a galaxy pair can¬ 
not affect the filament detection and the orientations of filaments 
and galaxy pairs are measured independently. 

The method and parameters for filament extraction are ex¬ 
actly the same as in Tempel et al. (2014a). The assumed scale 
(radius) of the filaments is roughly 0.7 Mpc'. Because the sur¬ 
vey is flux-limited, the sample is very diluted farther away, hence 
we are able to detect filaments of this scale only up to the redshift 
0.15 (=s 640 Mpc). To reduce the survey edge effects, we use only 
filaments that are farther than 10 Mpc from the survey boundary. 
In addition, to restrict ourselves with reliable filaments we dis¬ 
miss all filaments shorter than 10 Mpc. 


' Most of the galaxies are close to the filament axis and 90% of galax¬ 
ies inside 0.7 Mpc are closer than 0.5 Mpc to the filament axis. 
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Fig. 2. Upper panel', probability distribution function (blue solid 
line) of the projected (in the plane of the sky) angles between 
galaxy pairs and their host filaments. The KS-test value that the 
sample is drawn from a uniform distribution is 10“*. The hlled 
area shows the 95% confidence region for a randomised distribu¬ 
tion of 3012 pairs. Lower panel: the same as in the upper panel 
for two equal-size subsamples: close pairs (lisep < 0.3 Mpc; red 
dashed line) and loose pairs (^fsep > 0.3 Mpc; black solid line). 

2.4. Measuring the alignment signal 

As explained above, the determination of the orientation of fil¬ 
aments and the orientation of galaxy pairs are independent by 
construction and we do not have to worry about an intrinsic 
alignment caused by the method. However, the orientation of 
galaxy pairs and filaments are affected by the redshift-space 
distortions. To avoid any systematic effects introduced by the 
redshift-based distances, we measure the alignment in the plane 
of the sky. For that we project galaxy pairs and filaments to the 
plane of the sky and measure the angle between filament orien¬ 
tation and pair orientation. The measured angle is restricted to 
be in the range of 0 to 90 degrees. We calculate the alignment 
signal using a kernel density estimation.^ 

To estimate conhdence regions of the measured alignment 
signal, we generate 10000 Monte Carlo samples with the same 
number of pairs as in the original sample and conduct the same 
alignment analysis. In the figures below, we show the 95% con¬ 
fidence regions for this randomised signal. The Kolmogorov- 
Smirnov test (KS-test) probabilities that the angle distributions 
are drawn from a uniform distribution are also shown. To com¬ 
pare the two measured alignments, we use the KS-test to cal¬ 
culate the probability that the distributions are drawn from the 
same distribution. The KS-test values are indicated in the fig¬ 
ures. 

3. Results 

Figure 2 (upper panel) shows the probability distribution of the 
angles between galaxy pair orientations and their host hlament 
orientations. As explained above, the angles are measured in the 
plane of the sky in order to avoid possible effects of redshift- 
space distortions. It is seen that galaxy pairs are strongly aligned 
with galactic hlaments: the number of galaxy pairs aligned with 

^ In Appendix A of Tempel et al. (2014b) we show that the kernel 
density estimation is better than a simple histogram and is more repre¬ 
sentative of the underlying probability distribution. 



x(Mpc) (along filament) 

Fig. 3. The distribution of stacked galaxy pairs in the plane of the 
sky. Galaxy pairs are stacked based on the orientation of the host 
filaments (represented by the .r-axis); the centre point of each 
pair matches with the origin of the coordinate axes. The den¬ 
sities are shown in units of spherically averaged mean number 
densities. The black solid line shows the demarcation of close 
and loose pairs. 


the filament spine within +10 degrees is about 20% higher than 
would be in the case of a random distribution, and the pair num¬ 
ber decreases nearly monotonically when moving towards per¬ 
pendicularity with the filament spine. The KS-test probability 
that the alignment angles are drawn from a uniform distribution 
is Pks = 10“^ (~ 6.5cr signal). 

In the lower panel of Fig. 2 we have split the pair sample into 
two equal-size subsamples according to the projected distance 
between the galaxies of each pair. The demarcation between 
close and loose pairs (0.3 Mpc) is indicated in Fig. 1. Figure 2 
shows that galaxy pair alignment with galactic filaments is much 
stronger for loose pairs than for close pairs, while the alignment 
signihcance remains above 95% (2cr) in both cases. The KS- 
test probability that the close and loose pair samples are drawn 
from the same underlying distribution is 0.019, hence, the close 
and loose samples should be intrinsically different. However, the 
close-pair sample is contaminated by loose pairs oriented close 
to the line of sight, raising the measured amount of non-aligned 
pairs. 

We have checked also for the effect of sample splitting on 
the basis of various other measurable quantities: galaxy lumi¬ 
nosity, colour, magnitude difference, distance, morphology, dis¬ 
tance from filament axis, etc. In some cases, a small impact on 
the alignment signal was noticeable, but statistically insignih- 
cant. In principle, the statistical signihcance for the alignment 
signal can be increased either by increasing the sample size or 
by measuring the alignment in the real space rather than in pro¬ 
jection. The latter requires a proper treatment of redshift-space 
distortions, which is not a straightforward procedure and is be¬ 
yond the scope of the present analysis. 

To visualise the galaxy pair alignment in hlaments, we show 
a stacked distribution of the galaxy pairs in the plane of the sky 
in Fig. 3. For the stacking, the centre point of pairs was hxed 
at the origin of the coordinate axes and all hlaments were ori¬ 
ented along the x-axis. The hgure shows that for the loose pairs, 
the alignment occurs uniformly across the whole range of sep¬ 
arations. Note that the apparent density maximum in this hgure 
around [x = 0.35, y = 0.2] is actually statistically insignihcant - 
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Fig. 4. Upper panel: the alignment of close- and loose galaxy 
pairs with the host filaments in the Millennium simulation; i de¬ 
notes the real-space angle between a filament and a galaxy pair. 
Lower panel: the angle distribution for close and loose pairs in 
the SDSS data (solid and dashed lines; same as in Fig. 2) com¬ 
pared to the sky projection of a model distribution, in which 
the real-space alignment between pairs and filaments of the 
Millennium simulation (shown in the upper panel) is ascribed to 
the actual SDSS filaments. The modelled distributions are shown 
as filled regions corresponding to the 95% confidence level. 

it corresponds to an enhancement of the random distribution by 
just a few galaxies at the given projected distance. 

While the detected 2D alignment indicates that also in the 3D 
space, a significant fraction of galaxy pairs has to be aligned with 
filaments, a quantification and a more detailed interpretation of 
the measured signal is not straightforward due to projection ef¬ 
fects and redshift distortions. Instead, we can check whether the 
measured signal corresponds to the predictions of the ACDM 
formalism. To this end, using data from the Millennium simu¬ 
lation (as described above), we have studied the pair alignment 
with host filaments in the real space, splitting the pair sample 
into close and loose pairs according to the real-space distance^. 
The upper panel of Fig. 4 shows the resultant real-space align¬ 
ment with the filaments in the simulation data. As in the case 
of the SDSS data, both subsamples show significant alignment, 
with close pairs being somewhat less aligned than the loose ones. 
For a more direct comparison with the SDSS pairs, we have ap¬ 
plied the alignment of the Millennium pairs to the SDSS fila¬ 
ments. In 2D projection, the resultant angle distribution is shown 
in the lower panel of Fig. 4 as filled regions. For loose pairs, 
the pair alignment of the simulation matches well with the ob¬ 
served one. For close pairs, the distributions nearly match within 
the 95% confidence level, except for angles beyond about 70 de¬ 
grees, where the observed pairs seem to exhibit a population per¬ 
pendicular to the host filaments. However, the mismatch can be 
a result of the sample variance and does not provide grounds for 
invalidating the simulation. 

4. Conclusions and implications 

We used the SDSS Legacy survey to analyse the alignment be¬ 
tween galactic filaments and galaxy pairs. The filamentary pat- 

^ The use of real-space distances guarantees that the dependence 
on pair separation is intrinsic and not affected by projection effects. 
However, when splitting the sample based on the projected distance, 
the measured alignment signals remain roughly the same. 


tern in galaxy distribution was detected using the Bisous model 
and a catalogue of galactic filaments was constructed. To over¬ 
come possible contamination by galaxy pair orientations, the 
galaxy pairs were treated as single points during the filament 
extraction. The alignment is measured in the plane of the sky to 
avoid redshift-space distortions. 

The data show a clear tendency of galaxy pairs to align with 
their host filaments. The alignment is stronger for loose pairs 
and weaker for close pairs. The detected alignment is statisti¬ 
cally highly significant, reaching up to 6.5cr. Quantitatively, the 
alignment is similar to the alignment measured for galaxy pairs 
and filaments extracted from the Millennium simulation, thus 
being in concordance with the ACDM paradigm. 

Recently, the alignment of satellite galaxies and the cosmic 
web structures have been analysed by Lee & Choi (2015) and 
Tempel et al. (2015). Both studies conclude that the distribution 
of satellites around primary galaxies is elongated in the direction 
of galactic filaments. 

A simulation-based test carried out by Tempel et al. (2014c) 
showed that the given filament detection method also recovers 
the underlying velocity field. We can thus deduce that the align¬ 
ment of galaxy pairs and satellite systems with their host fila¬ 
ments is a direct result of the anisotropic velocity field, leading 
towards a preferred direction of satellite accretion and gravita¬ 
tional collapse. 
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